Giant electrocaloric effect around T c 
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We use molecular dynamics with a first-principles-based shell model potential to study the elec- 
trocaloric effect (ECE) in lithium niobate, LiNbC>3 , and find a giant electrocaloric effect along a 
line passing through the ferroelectric transition. With applied electric field, a line of maximum 
ECE passes through the zero field ferroelectric transition, continuing along a Widom line at high 
temperatures with increasing field, and along the instability that leads to homogeneous ferroelectric 
switching below T c with an applied field antiparallel to the spontaneous polarization. This line is 
defined as the minimum in the inverse capacitance under applied electric field. We investigate the 
effects of pressure, temperature and applied electric field on the ECE. The behavior we observe in 
LiNbC>3 should generally apply to ferroelectrics; we therefore suggest that the operating tempera- 
ture for refrigeration and energy scavenging applications should be above the ferroelectric transition 
region to obtain large electrocaloric response. We find a relationship among T c , the Widom line 
and homogeneous switching that should be universal among ferroelectrics, relaxors, multiferroics, 
and the same behavior should be found under applied magnetic fields in ferromagnets. 



PACS numbers: 77.70.+a,77.80.-e,34.20.Cf,83.10.Rs 

The electrocaloric effect (ECE) has the potential for 
efficient solid-state refrigeration pQ, particularly for chip- 
level applications, but fundamental questions about the 
nature of the ECE have plagued the field for at least 50 
years [21 [3]. Here we show how the ECE can be large 
well above T c , in spite of the apparent thermodynamic 
argument that it vanishes above T c . The ECE is the 
change in entropy with electric field, and thus can be 
used to pump heat or scavenge energy. Experiments 
show temperature changes of up to 12 K on applica- 
tion of reasonable electric fields of 480 kV/cm in PZT 
[H [5] . Simulations with an effective Hamiltonian of the 
ECE show good agreement with experiments, [3 [7] al- 
though there are questions about the observed behavior 
above T c , discussed below. We study the behavior of 
LiNb03 as a simple, well-studied, uniaxial ferroelectric 
[SI [3] for which we have a well-tested first-principles based 
interatomic potential 10- 13 , with which we previously 
studied the converse effect, pyroelectricity, under zero 
applied field[13]. Recently the ECE in Ba .5Sro.5Ti0 3 
was simulated using an effective Hamiltonian at constant 
volume. [14] Here we use an interatomic potential which 
includes all degrees of freedom to perform simulations at 
constant pressure and temperature, and develop a differ- 
ent picture of the ECE from Ref. El as discussed below. 
We carefully study the ECE around T c , and find that 
ferroelectric switching below T c is closely related to the 
polarization response above T c under applied fields. 

The thermodynamic inverse of the pyroelectric effect, 
the electrocaloric effect is the reversible temperature 
change of a polarizable material under the application 
and/or removal of an electric field. Previous studies have 
focused on temperatures below the transition tempera- 
ture from the ferroelectric state to the paraelectric state. 
These studies have also kept their focus on applied fields 



in the same direction as the spontaneous polarization of 
the ferroelectric state. 

We performed molecular dynamics (MD) simulations 
on supercells of 8x8x8 unit cells, with a total of 5210 
atoms using the DLPOLY package [15] with a Nose- 
Hoover thermostat and barostat in the NaT ensemble. 
The shell model potential was fit to first-principles den- 
sity functional theory (LDA) energies, forces, stress, ef- 
fective charges, dielectric constants, phonons frequencies 
and eigenvectors and has proven predictive and robust. 
At lower temperatures (relative to T c ) very long runs 
were needed to reach equilibrium; e.g. at 800 K, 22.5 ns 
was required to equilibrate the shell and atomic positions. 
Once equilibrium was reached, the stable configuration 
was used as a starting point for runs at similar temper- 
atures. For all simulations with the electric field applied 
in the direction opposite to the polarization, the simula- 
tions were cither 44 or 110 ps long with a timestep of 0.2 
and 0.5 fs, respectively. At 4 or 10 ps intervals the electric 
field was increased by 5 MV/m. For simulations with the 
field applied parallel to the spontaneous polarization, we 
ran a simulation at a fixed applied field and temperature 
for 100-200 ps with a timestep of 0.5 fs. All averages were 
taken after the equilibration period, which was typically 
25-50 ps. Simulations with applied pressure ran for 100 
ps with a timestep of 0.5 fs. We employ the convention 
defining fields applied in the same direction as the spon- 
taneous polarization as positive or parallel and the fields 
applied in the direction opposite that of the spontaneous 
polarization as negative or antiparallel. The results of po- 
larization versus temperature and applied field are shown 
in Fig. 1. We find the zero-field ferroelectric to paraelec- 
tric transition to be 1506K, which slightly overestimates 
the literature value of 1480K[S]. 

To investigate the electrocaloric effect we used the 
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FIG. 1. (a) Average polarization of LiNbOa with applied field 
along the same direction as the spontaneous polarization and 
(b) derivative of polarization with respect to temperature. 
The Widom line T* above T c is traced out by the inflections 
in P or minima in dP/dT. 
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where P is the macroscopic polarization, E is the macro- 
scopic electric field, C is the heat capacity, V is the vol- 
ume, T is the absolute temperature, and p is the pressure, 
which when integrated gives the thermodynamic identity 
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to calculate the change in the temperature. The maxi- 
mum AT is found at progressively higher temperatures 
as the strength of the field is increased (Fig. 2). We also 
performed constant energy (microcanonical, NVE) sim- 
ulations with applied fields, and got agreement within 
error to the results in Fig. 2, obtained from equation 3, 
considering the differences between constant V and con- 
stant p simulations. 

For small negative or antiparallel fields below the fer- 
roelectric transition we expect to find the same values of 
the ECE as for small positive fields since (j^)\e is con- 
tinuous around E = 0. When we apply a negative field 
larger than the coercive field E c at a given temperature 
Ths , the polarization switches homogeneously, and the 
direction of the polarization of the crystal becomes par- 
allel with the direction of the applied field (see page 79- 
80 in Ref. [TH] for a discussion of this instability) . We call 
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FIG. 2. Calculated change in temperature as a function of 
temperature for various applied electric fields. AT above 
T c increases with applied field since the polarization P is 
higher. Note that points shown here are calculated and curves 
are fit to the data. 



this temperature for a given electric field the homogenous 
switching temperature, Th s ■ Normally switching occurs 
through growth of preexisting domains with a more pre- 
ferred polarization direction, or nucleation and growth of 
new domains. Homogeneous switching, with no preexist- 
ing domains, is the limit of metastability of the energeti- 
cally unfavored polarization direction in an applied field. 
We recently studied homogeneous switching in detail in 
PbTiC>3 [T7] . This switching occurs, depending on the 
strength of the applied negative field, on the time scale 
of our simulations between 4 and 10 ps. 
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FIG. 3. Average polarization for LiNbOa with applied elec- 
tric field in the opposite direction to spontaneous polariza- 
tion. The negative ordinate represents the direction of the 
polarization after switching to be parallel to the applied field. 

We performed constant energy (NVE) simulations 
to study the temperature change during polarization 
switching. As expected, the temperature goes up sharply 
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with switching: since the potential energy drops, the ki- 
netic energy, and thus T, increases. By varying the ap- 
plied field frequency, however, it should be possible to get 
the polarization switching to go 180° 's out of phase and 
make a refrigeration cycle based on polarization switch- 
ing. Such behavior may be an explanation of the depen- 
dence of ECE on applied field frequency seen in experi- 
ments pQ. 

In order to understand the line of maximum ECE, 
we consider the thermodynamics of ferroelectrics in an 
applied electric field. At zero field the ferroelectric 
transition, Ti , is first-order, and with applied field, 
Ti increases as = — 4^ [15] . With increasing ap- 
plied field, a tricritical point, T tr i , is reached where the 
transition becomes second-order, T s . In BaTi03 this 
first-order region between and T tr ; is 18K. 3 With a 
further increase in field a critical point T c is reached, and 
there is no phase transition above that temperature and 
applied field, E c . We were not able to resolve the re- 
gion between Ti and T c in our MD simulations at zero 
pressure. Although there is no phase transition above 
T c for higher fields, there is still an inflection point in P 
versus T and a line is traced out with increasing E that 
can be considered the Widom line for this problem [T§1 - 
[22] . We label the Widom line T* . In a perfect crystal 
there is an instability with increasing negative field at 
each temperature below Ti , so negative fields greater in 
magnitude than the coercive field will switch the polar- 
ization. This traces the line of homogeneous switching, 
Ths • This instability is found clearly in first-principles 
computations at constant electric field and zero temper- 
ature in PbTiC>3 (see Fig. 6c in Ref. |2"3"]) . The impor- 
tance of the inverse capacitance C _1 oc e" 1 , which is the 
second derivative of the internal energy with electric dis- 
placement field D, where e is the dielectric constant as 
discussed in Ref. [24] We propose that a line of minima 
in the inverse capacitance C _1 joins T^ , Ti , T tr j , T c , 
and T* . Below T c , C _1 = along this line, whereas it 
is a minimum above T c defining T* . These temperatures 
are sketched in the inset to figure 5 below. 

As J.F. Scott discusses, the issue of ECE above T c has 
been a argument spanning many decades over whether 
the ECE is expected above or below T x , and whether 
it has to be related to a ferroelectric phase transition [3] • 
The confusion stems from misconstruing the Maxwell re- 
lation (Eq. 1). Above Ti , P = at zero field, implying 
zero ECE, but at larger fields ? ^ so that there is 
an ECE in any material, not just a ferroelectric. Fur- 
thermore there is no requirement of an underlying fer- 
roelectric phase transition. The fact that the dielectric 
constant e _1 oc C _1 is high above Ti in a ferroelectric, 
however, makes the application of ferroelectrics above 
Ti optimal, rather than far below Ti as seems to be stan- 
dard practice now. The susceptibility is maximal along 
the Widom line, so that is where we expect to see the 
greatest ECE, as we find. The maximum in AT occurs 



along the line defined by minC -1 , and our results clarify 
the issue of ECE above Ti . The order parameter Q for 
polarization correlations Q — lim. r _ i . 00 < P(r)P(0) > \e, 
and the maxima in the ECE should be where dEdT is 
maximal. 

In contrast, Ponomareva and Lisenkov state that "AT 
peaks near the transition temperatures which is the con- 
sequence of the larger configurational disorder associated 
with these points." [2] . To the contrary, we find that AT 
continues to increase as the electric field increases, well 
above the zero field Ti . So we find that operation above 
Ti would be advantageous in ECE applications. The is- 
sue of disorder will be discussed further below. 

We have also simulated LiNbC>3 as a function of applied 
pressure. As pressure increases, Ti strongly decreases, so 
that it is 750K at 40 GPa. As pressure increases and 
Ti decreases, we find consistency with the picture above, 
with the ECE maximal along the Widom line. Unexpect- 
edly, we found that pressure increases the phase transi- 
tion region between Ti and T tr i under applied field, so 
that it is more easily discernible (Fig. 4). We also find 
a structural phase transition between 60 and 70 GPa, 
broadly consistent with the transition seen experimen- 
tally at about 40 GPa[Hl[27]. We did not attempt to 
reverse this transition for this study of the ECE, but it is 
encouraging that our potential does reproduce the high 
pressure transition. 




FIG. 4. Effect of pressure and electric field on polarization in 
LiNbOa . Note the strong decrease in the transition temper- 
atures with pressure, and the increase in the phase transition 
region which appears to be first order up to at least 25 MV/m 
and about 100K at 20 GPa. At zero pressure the phase tran- 
sition disappears at relatively small fields. 

Our results are consistent with the above discussion 
and agree well with recent experiments. For example, 
Bai et al. found that in BaTiOa films the maximum ECE 
occurs above T c and increases with applied field [35] [23] . 
Relaxors and relaxor ferroelectrics are expected to show 
similar behavior to what we see in LiNbOs (e.g. Fig. 3 in 
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FIG. 5. Effects of pressure and electric field on the transi- 
tions in LiNbC>3 . Ti , T c , T* along the Widom line, and 
Ths , the homogeneous switching temperature (see text) fall 
along a smooth curve. The right hand figure shows the line 
of minimum inverse capacitance, and the inset shows the re- 
lationship between the various temperatures along these min- 
ima schematically, as discussed in the text. We could not 
constrain the location of Ttri and T c is our simulations at zero 
pressure due to the small temperature region between them 
and Ti . is presented by the dashed red segment, the first- 
order transition region between Ti and T tr i is the solid black 
segment, T s is the dashed blue segment between T tr i and T c , 
and T* is the dotted violet high T segment. 



Ref.[H]). 

One might find it surprising at first that an apparently 
kinetic phenomenon, ferroelectric switching, is related to 
the thermodynamics at Ti , and moreover, to the behav- 
ior high above the zero field T c where there is no phase 
transition at all. We can understand this behavior from 
an atomistic perspective by considering a simple classical 
double well representing the underlying potential surface 
(Fig. 6). The key is that the same double well poten- 
tial surface underlies ferroelectric switching, Ti , and the 
maximum in e oc C: In zero applied field, Ti is given by 
coupled local double well potentials. In the displacive 
limit Ti is the temperature where there is enough kinetic 
energy to traverse the hump in the underlying effective 
double wells. In an applied field the underlying poten- 
tials are asymmetric, and it takes a higher T to go over 
the hump. There is no true phase transition above the 
critical field, but the underlying potential surface is still 
asymmetric due to the field, and the susceptibility, which 
diverges at a true phase transition, instead peaks at the 
temperature where local modes start to hop over the 
hump. This is the Widom line. Below Ti the effective 
potentials are asymmetric due to the coupling between 
local modes. In other words the coupling can be thought 
of as a local Weiss field from the rest of the crystal. When 
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FIG. 6. The connection between the transition temperatures 
and switching below Ti with an antiparallel applied field is 
shown in a simple schematic. These do not show the Landau 
free energy, but rather the potential energy versus local polar- 
ization. In all cases an activation barrier exists for changing 
the local polarization. 



a field greater than the coercive field is applied, the po- 
larization switches. In terms of the double well picture in 
figure 6, the field makes a population inversion, where the 
high energy state is occupied. At low temperatures or be- 
low the coercive field for a given temperature the atoms 
will be stuck in the high energy state, for a time related 
to the hump size and the asymmetry of the double wells 
(the driving force). As temperature is raised for a given 
time scale there will be a point at which the local modes 
all cross the hump and go into the lower energy state. So 
we see the origin of the connection between switching, 
Ti , and the Widom line. In the order-disorder limit, 
with a huge hump that cannot be broached, the dynam- 
ics are governed by nucleation and growth, but we still 
expect to see a relationship between switching, Ti , and 
the Widom line since the field dependence of the driving 
force for the transition and for switching is still governed 
by the E ■ P term in the free energy. Finally, pressure 
reduces Ti by decreasing the ferroelectric well-depths. 

This behavior for applied fields carries over as well to 
magnetism, and surprisingly does not seem to have been 
remarked upon in the literature. In experiments for mag- 
netization of pure iron under applied magnetic fields, the 
Widom line is clearly seen about T c , and the Widom 
line comes in with a finite slope at T c [3U] . This suggests 
that magnetic switching should also follow this line below 
T c , otherwise the Widom line would come in with zero 
slope. So homogeneous switching, T c , and the Widom 
line in ferromagnets is related in the same way to ap- 
plied magnetic fields, as ferroelectrics behave in applied 
electric fields. 

Some have attributed the ECE to polarization disor- 
dering giving a resulting change in entropy |14) , but this 
simple view of statistical mechanics can be misleading. 
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For example, an isolated supercooled liquid will crystal- 
lize with a concomitant increase in entropy, as required 
by the second law of thermodynamics, even though a 
crystal is more ordered than a liquid. Pressure decreases 
T c , and thus increases the ECE, not because it directly 
causes greater disorder (an entropic effect), but because 
of an cnthalpic, volumetric effect: the paraelectric phase 
is denser than the ferroelectric phase, and thus favored by 
pressure. The Maxwell relation, equation 1, says every- 
thing necessary about the relationship of the ECE to the 
pyroelectric effect. Furthermore, this is an exact expres- 
sion, and thus the "direct" and "indirect" approaches to 
the ECE must agree within the error of experiments or 
simulations, in contrast to their discussion. It is also im- 
portant to realize that the A.S tr , the entropy change for 
the ferroelectric to paraelectric first-order phase transi- 
tion is not the pertinent quantity for the ECE, but rather 
(dS/dE)\T, which is not the same thing. Large AS tr can 
give rise to a large ECE near Ti , but only in the first- 
order transition region. We find larger ECE at higher 
fields and temperatures beyond T tr i ■ 

Our MD simulations with first-principles potentials 
on LiNb03 give some universal conclusions about the 
electrocaloric effect. We expect the electrocaloric ef- 
fect (ECE) to be maximal above Ti , suggesting ma- 
terials with Ti below room temperature (or the oper- 
ating temperature) would be optimal. We find a rela- 
tionship between Ti , the Widom line, and homogeneous 
switching, and suggest this may be universal among fer- 
roelectrics, relaxors, multiferroics and ferromagnets. The 
ECE should be large in any insulator with a large, tem- 
perature dependent, dielectric susceptibility, not just fer- 
roelectrics. We have shown the importance of the Widom 
line and the nature of the ECE above Ti , the rela- 
tionship between homogeneous switching and behavior 
above Ti in an applied electric field, the continuity of 
the Widom line with the switching behavior below Ti , 
the realization that ECE does not require a ferroelec- 
tric; the relationship of switching and the Widom line in 
ferromagnets as well as ferroelectrics, clarification of con- 
ditions to maximize the ECE in materials, the behavior 
of the ECE under pressure, and finally the application of 
an accurate first-principles based shell-model to behavior 
in an applied electric field. 
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